The Charmonium (Bottomonium) binding ate finite T 
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The charmonium (bottomonim) binding at finite temperature is studied with static potentials 
extracted from the lattice QCD data of Kaczmarek et al. The bottomonium spectrum is also 
studied. This is relevant for Hard Probes in Heavy Ion Collisions. 
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I. INTRODUCTION 

This talk is motivated by the seminal hard cc probe 
paper of Matsui and Satz [1[ by the lattice QCD data on 
finite temperature static potentials [2, y, 0, 0, H S S, @| 
and bindin g of charmonia [13, \lM > ^^id by the theoreti- 
cal studies [l^, [l^ [ij, [l^ on finite temperature lattice 
potentials. 

With a finite T quark potential, and modern quark 
model techniques one might, approximately, 

- study chiral symmetry breaking, quark mass genera- 
tion, a finite T notice that if one maintains a confining 
potential, chiral symmetry is always broken, whatever 
the T, 

- compute the spectrum of any hadron at finite T not 
only the J/tp but also light mesons, baryons, etc 

- compute the interaction of any hadron-hadron at fi- 
nite T using cluster methods like the Resonating Group 
Method. 

Here just study the charmonium and bottomonium as 
prototypes to study finite T quark potentials. 

The charmonium is a good starting point because 
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Thus it is reasonable to neglect in the bound state equa- 
tion, spontaneous chiral symmetry breaking, relativistic 
effects, coupled channels, and temperature effects (other 
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FIG. 1: Lattice QCD data for the internal energy U\, thanks 
to Olaf Kaczmarek et al. The solid line represents the T = 
potential . 
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FIG. 2: Lattice QCD data for the free energy F\, thanks to 
Olaf Kaczmarek et al. The solid line represents the T = 
potential . 



than the potential dependence on the temperature) . We 
can simply solve the Schrdinger equation with static lat- 
tice QCD finite T potentials 

Notice that a Coulomb potential is sufficient to bind 
an infinite number of charmonia, and therefore the loss of 
the linear confinement is not sufficient to melt the char- 
monia states. A detailed fit of the finite temperature 
potentials is needed to determine the melting of charmo- 
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FIG. 3: Boundstates with the Free Energy potentials. In the top left we show the energy spectrum of charmonium in the 
enveloping potential, and in the top right we show the wavefunctions for the ground state. In the bottom left we show the 
energy spectrum of bottomonium in the enveloping potential, and in the bottom right we show the wavefunctions for the ground 
state. 



II. THE FINITE TEMPERATURE STATIC 
QUARK POTENTIALS 



We assume that the static quark-antiquark potential V 
is due to the length r of the flux tube, related the volume 
Vvoi of excited QCD vacuum confined into the flux tube. 
Our first step consists in assessing V ^ the static quark- 



c 


a 


a' 


A 


O" 


-0.0270 


0.2857 


0.1333 


10.00 


1.1621 


TjT, 


d 


A 


A 


w 



1.05 


0.98 


0.357 


2.1 


0.10 


1.20 


0.70 


0.340 


2.7 


0.05 


1.50 


0.50 


0.537 


2.7 


0.10 


3.00 


-0.045 


0.301 


6.0 


0.05 


6.00 


-0.68 


0.239 


11.0 


0.10 


12.0 


-1.55 


0.116 


25.0 


0.10 



antiquark potential, 



dV 



-adr , 



c a a A a 


-0.2117 0.3490 0.1733 7.793 4.475 


T/Tc d A X A' X' 


w 



1.13 


3.94 


6.99 


2.64 








0.336 


1.18 


2.35 


4.58 


2.27 








0.847 


1.27 


2.64 


5.38 


2.87 








0.461 


1.40 


2.23 


5.57 


2.91 











1.64 


1.79 


5.60 


3.65 











1.95 


1.51 


5.57 


4.15 








0.117 


2.61 


1.21 


6.12 


5.44 











4.50 


-40.7 


12.6 


5.44 


1.00 


0.018 





7.00 


-74.0 


15.9 


5.44 


0.49 


0.021 






(2) 



TABLE I: Fitting parameters of the free energy Fiy(r) 



TABLE 11: Fitting parameters of the internal energy Uirir) 
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FIG. 4: Boundstates with the Internal Energy potentials. In the top left we show the energy spectrum of charmonium in the 
enveloping potential, and in the top right we show the wavefunctions for the ground state. In the bottom left we show the energy 
spectrum of bottomonium in the enveloping potential, and in the bottom right we show the wavefunctions for the ground state. 



where here cr is a force, generaUzing the string tension. 
In lattice QCD, with the Polyakov Loop, one computes 
Fi, the Free Energy 



dFi ^ -a dr - S dT 



(3) 



which equals the potential for isothermal transforma- 
tions. To compute Ui, the internal energy, one needs 
to compute the entropy S with the total action as well 



dUi = —crdr + TdS, 
= dFi - d{TS) , 



(4) 
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which equals the potential for adiabatic transformations. 
Naturally close to the phase transition temperature Tj,, 
the transformations are nearly isothermal and, and the 
potential V is close to the internal energy Ui, while far 
from Tc the transformations are nearly adiabatic poten- 
tial is close to the free energy Fi . The interpolating rela- 
tion for V between Fi and Ui has been derived by Wong. 

Lattice QCD provides potentials and energies, both 
with quenched (pure gauge) and dynamical (with 
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TABLE III: Goundstate solutions of the charmonium at dif- 
ferent temperatures, for the free energy i^iy(r) . 



TABLE IV; Groundstate solutions of the bottomonium at 
different temperatures for the free energy Fixir) ■ 



fcrmions) for quarks in confined system such as mesons, 
diquarks, baryons, tetraquarks, pentaquarks and hy- 
brids. Lattice also computes spectral distribution func- 
tions, to access directly, say, the bottomonium masses. 
However the limitations of the lattice QCD potentials 
are, 

- static potentials only have been computed, 

- a constant shift of the potential is undetermined 

- few spin dependent potentials are so far computed, 

- lattices have a relatively small volume. 
Nevertheless we can extract the potentials from lattice 
data and use them to compute the spectra with the 
Schrodinger equation. 



distance par of the energy. 

We fit the enveloping function V{r) of the free energies 
Firp{r) with a constant shift c, two Coulomb potentials, 
one screened for the short distance and another for the 
long distance part of the potential, and a linear potential. 



V{r) = c + 



-a + a'e-^'' 



(5) 



The parameters are shown in Table HI 

We fit the long distance part of the free energies Fixir) 
with a saturation function Mxir), including a constant 
shift d and an exponential decay, 
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III. FITTING BETWEEN THE T=0 ENVELOPE 
AND THE FINITE T SATURATION 



In the plots for the free energy _Fiy and the internal 
energy Uit for the different temperatures T as in Figs. [T] 
and m it is clear that each set of energies is bound from 
above by a common enveloping function. 

The enveloping functions are respectively the T = 
energies, Uiq and Fig. Notice that the T = free (inter- 
nal) energy also coincides with the small distance part of 
the free (internal) energy functions. 

In what concerns the large distance par of the free 
(internal) energy, for T < Tc the string tension decreases, 
while for T > Tc the string tension vanishes and the free 
(internal) energy saturates. 

Thus we fit the free (internal) energies in three steps: 

- fist we fit the enveloping function 

- then we fit the long distance saturation 

- finally we match the long distance part to the short 
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The parameters are shown in Table IH 

The interpolation between the short distance Fo{r) and 
the long distance Mxir) is performed with the coupled 
channel method, where we choose the lowest eigenvalue 
of the matrix 



V{r) wt 
wt* Mrir) 



(7) 



where the interpolating parameter wt is shown in Table 
III 
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TABLE V: Charmonium boundstates computed with the in- 
ternal energy UiT{r) 



TABLE VI: Bottomonium boundstates computed with the 
internal energy t/iT(r) 



Finally the free energy is, 



V{r) + Mrir) - J[V{r) ~ Mrir)]^ + |wp 
Frir) = ^— ^ . 

(8) 
The fit is achieved with 5 constant parameters for V{r), 
and 4 parameters per temperature for the long distance 
part Mxir) and for the matching parameter w. 

In what concerns the internal energy Uix{r), a fit with 
similar functions but with different parameter is per- 
formed, except for the higher temperatures that need an 
extra pair of parameters in the saturation function. 
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The parameters of the internal energy are listed in Table 

im 



IV. SOLVING THE SCHRODINGER EQUATION 



To solve the Schrodinger equation 



-—V^ipir) + V(rU = Eil;(r) 
2m 
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TABLE VII: Bottomonium boundstates computed with the 
internal energy C/it('') (continued) 



with a potential depending only in r, we can separate the 
eigenf unctions as 



^{v) = ^YUe,v), 



(11) 



and we get the following equation for the radial compo- 
nent u{r) {Ti — 1 ) 



I (Pur^l , l{l + l) 



2m dr^ 
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Unl + V{r)Unl == EnlUnl ■ (12) 



To solve the radial equation we discretize it, with the 
finite difference substitution 
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We also impose the Dirichlet boundary conditions uo = 
um — 0, consistent with the radial equation. So we get 
the linear system 



H.jUj 
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(14) 



where Hij is tridiagonal. We solve the linear system 
for the lowest eigenvalues, by using the inverse iteration 
method. 

The results depend on the mass of the heavy quarks, 
affecting the kinetic energy and on the temperature, af- 
fecting the potential. The charm mass and the bottom 
mas, in different potential models, ranges respectively 
from 1300 MeV to 1752 MeV and from 4750 to 5100 
MeV. 

Since the lattice potentials have a constant energy shift 
undefined, we don't show the total energy of the systems, 
only the binding energy, defined with the difference be- 
tween the potential at infinity and the boundstate en- 
ergy. For instance for T = 0, where confinement occurs, 
the binding energy is — oo. 



V. RESULTS 

We now show the results of the boundstate equations. 

A fisrt study of the melting of the boundstates can 
be achieved with the enveloping potentials, for the in- 
ternal energy and for the free energy. Comparing the 
energy levels of the enveloping potentials with the satu- 
ration energies at the different temperatures, we estimate 
the melting temperatures. Essentially, we can estimate 
melting to occur when the energy level is close to the 
saturation energy. 

Thus we first solve the Schrodinger equation with the 
V{r) enveloping potential for the charmonium and for 
the bottomonium and compare the energy levels with 
the potentials obtained in finite temperature. This is 
depicted in Figs. [3] and ID It occurs that only few of 
the tempetarures T > Tc will provide binding for the 
charmonia, while the bottomonia survives up to higher 
temperatures. It also appears that the internal energy 



may provide binding up to higher temperatures than the 
free energy. 

We then study binding for the different potentials 
extacted fom lattice QCD simulations of Kaczmarek 
05 S IJj S 0- with the free energy (quenched or dy- 
namical) the J/'ip and rjc melt at T/Tc ~ 1.17 and the 
groundstate bottomonium melts at T/Tc — 1.8. With 
the internal energy (quenched) the groundstate charmo- 
nium melts at T /T^ ~ 1.6, while the p-wave excitations, 
say the Xcj nielt at T/Tc — 1.1, while the bottomonium 
groundstate melts at T/T^ ~ 3.2, the p-wave Xb melts at 
T/Tc =i 1.7, the d-wave melts at T/T^ ~ 1.2. Examples 
of wavefunctions are also show in Figs. [3]and|4l The de- 
tails of the binding are shown in Tables ITTTl and ITVl for the 
free energy and in Tables fV l I VIII and I VIII for the internal 
energy. 



VI. CONCLUSION 



study reproduces the charmonium melting temperatures 
in the literature. 

In particular the J/ij] melts above Tc (at 1.6 Tc ac- 
cording to lattice results) but the excited charmonium 
Xc or -0* are probes that melt just above Tc. New in our 
results are the detailed studies of the bottomonium, rele- 
vant for the LHC, with a melting as high as T/Tc — 3.2. 

Quantitative Puzzles remain in understanding the the 
mass shifts > 300 MeV , unseen in the experimental data, 
and in the importance of the the spin dependent poten- 
tials, missing in the lattice QCD data. 

Possible future efforts may be to, 

- calibrate the best we can the T-dependent quark- 
antiquark potentials, 

- apply the potentials to chiral symmery break- 
ing/restoration, 

- apply chiral restoration to the light quark mesons p, w, 4> 
and TT, K. 



Wong [l^l studied how the static potential interpolates 
between the free energy and the internal energy at dif- 
ferent temperatures. Wong showed that at T ~ T^ the 
free energy approximates the static potential, while at 
T » Tc the static potential is closer to the internal 
energy. 

We study the binding of charmonium and bottomo- 
nium, both for the free energy and for the internal en- 
ergy, with a finite difference lattice of 200001 points. Our 
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